Kosterlitz-Thouless behavior in layered superconductors: 
the role of the vortex-core energy 
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In layered superconductors (SC) with small interlayer Josephson coupling vortex-antivortex phase 
fluctuations characteristic of quasi two-dimensional (2D) Kosterlitz-Thouless behavior are expected 
to be observable at some energy scale Td. While in the 2D case Td is uniquely identified by the KT 
temperature Tkt where the universal value of the superfluid density is reached, we show that in a 
generic anisotropic 3D system Td is controlled by the vortex-core energy, and can be significantly 
larger than the 2D scale Tkt- These results are discussed in relation to recent experiments in 
cuprates, which represent a typical experimental realization of layered anisotropic SC. 

PACS numbers: 74.20.-z, 64.60.Ak, 74.72.-h 



Since the pioneering work of Kosterlitz and Thouless[l] 
(KT) on the so-called KT transition in the two- 
dimensional (2D) XY model, much attention has been 
devoted to the effect of phase fluctuations in quasi 2D 
superfluid systems. Thin films are natural candidates 
for the observation of KT physics, as the occurence of 
the "universal" (i.e. sample-independent) jump of the 
superfluid density, measured in ^He superfluid films, or 
the non-linear I — V characteristic, observed in thin 
films of conventional SC. [2] Signatures of KT physics 
can be expected also in layered SC with weak inter- 
layer coupling. A remarkable example of systems be- 
longing to this class are underdoped samples of high- Tc 
SC. [3] Recently, various experiments ranging from finite- 
frequency conductivity, [4, 5] Nerst effect [6] and non- 
linear magnetization[7] have been interpreted as signa- 
tures of KT phase fluctuations. Nonetheless, any effect 
reminiscent of the universal jump of the superfluid den- 
sity at Tkt, which would be the most direct probe of KT 
physics, failed to be observed[8-12]. 

Until now, the 2D-3D crossover in anisotropic layered 
SC has been discussed mainly within the framework of 
the anisotropic 3D XY model[I3-17] 
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Here 



is the superconducting phase on two nearest- 



neighbor sites (i,j) of a coarse-grained lattice, on the 
same plane (Jij = Jab) or in neighboring planes (Jy = 
Jc). The energy scales Jab, Jc can be related to the mea- 
sured 3D superfluid density ps at T = as: 
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where Aa6,Ac represent the in-plane and out-of-plane 
penetration depth, respectively, m is the electron mass, a 



is the in-plane lattice spacing and d is a transverse length 
scale (i.e. the interplane distance) used to define the ef- 
fective 2D areal superfluid density p^'^ = dps- In a 2D 
system vortex fluctuations drive p1'^{T) to zero at Tkt 
given by (we put h = c = ks = I) 
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where the temperature dependence of pj. {T) includes 
also the effect of other excitations, like long- wavelength 
phase-fluctuations of the model (1) or BCS-like quasi- 
particles excitations. [3, 9, IS, 19] Within the anisotropic 
XY model (1) a flnite interlayer coupling Jc cuts off the 
logarithmic divergence of the in-plane vortex potential 
at scales ~ a/Y^,[13] where 77 = Jc/ Jab, so that the su- 
perfluid phase persists above Tkt, with Tc at most few 
percent larger than Tkt- [14-16] As far as the superfluid 
density is concerned, there is some theoretical[14] and 
numerical [17] evidence that even for moderate anisotropy 
the universal jump at Tkt is replaced by a rapid down- 
turn of Ps{T) at a temperature scale Td — Tkt- 

However, recent measurements of Ps{T) in strongly un- 
derdoped YBa2Cu306+:r (YBCO) samples[ll, 12] (with 
large 77 ^ lO^'* anisotropy [12, 20]), showed that no down- 
turn of Ps (T) is observed at the KT temperature defined 
by Eq. (2), but eventually at a scale Td « ?c.[ll] Anal- 
ogously, recent measurements of the phase-fluctuations 
diamagnetism in underdoped Bi2212[7] revealed that the 
phase correlation length ^ above Tc can be fltted with 
the typical KT law, provided that the effective KT tem- 
perature is few kelvin smaller than Tc- In both cases, by 
looking at the system from below or above Tc, it appears 
that the typical temperature scale where vortex fluctua- 
tions become relevant is always near Tc, regardless of the 
value (2) of the Tkt of the pure 2D case. 

In this paper we analyze the role played by the inter- 
layer coupling and the vortex-core energy at the crossover 
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from 2D KT to 3D superconducting behavior in layered 
SC. In particular, we focus on the behavior of the super- 
fluid density below T^. and of the correlation length above 
Tc- We carry out a rcnornialization group (RG) analy- 
sis using the mapping between the thermal metal-SC KT 
transition in 2D and the quantum metal-insulator tran- 
sition in the ID sine-Gordon model. [22] Indeed, a similar 
model has been studied in Ref. [21] to investigate the 
superfluid- insulator transition in optical lattices of ID 
boson chains, where the tunneling inter-chains amplitude 
plays the same role of the Josephson coupling in layered 
SC. We show that in the presence of a finite interlayer 
coupling the superfluid density looses its universal char- 
acter. The jump in Ps{T) at Tkt observed in the 2D case 
is replaced by a downturn curvature at a temperature 
which depends on the vortex-core energy /i. While in 
XY models, where n is fixed by the in-plane coupling 
Jab (see Eq. (6)), ~ Tkt, in the general case the ratio 
Tcl/Tkt increases as j^i/Jab increases. Analogously, by 
approaching the transition from above, the increasing of 
the phase-fluctuation correlation length is controlled by 
the scale Td instead of the Tkt of the pure 2D system. 
Based on these results, we argue that the various exper- 
imental data in cuprates concerning KT behavior can be 
reconciled if /i is larger than the typical XY value. 

Let us first recall briefly the basic features of the KT 
transition using the analogy with the quantum ID sine- 
Gordon model, [2, 22] defined as: 
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Following standard definitions, [22] 6 and dx4> represent 
two canonically conjugated variables for a ID chain of 
length L, with [9{x'), dx(l){x)] = i'KS{x' — x), K is the 
Luttinger-liquid (LL) parameter, Vs the velocity of ID 
fermions, and Qu is the strength of the sine-Gordon po- 
tential. For g„ = the action of the model (3) can be 
simplified by integrating (f) and rescaling r — > VsT, so that 



So = ^J dxdr [{d^ef 



(4) 



equivalent to the gradient expansion of the model (1), 
with T as the second spatial dimension. Note that in the 
Hamiltonian notation (3) the coefficient of the dual field 
is 1/K, while the rotational in-plane symmetry of the 
model (1) is recovered in the action (4). Besides the long- 
wavelength phase fiuctuations present in Eq. (4), vortex 
configurations are possible, which require §V0 = ±27r 
over a closed loop. Since (j) is the dual field of the phase 
6, a 2iT kink in the field 9 is generated by the operator 
g«20 j22] i.e. by the sine-Gordon potential in the Hamilto- 
nian (3). The correspondence between the quantum ID 
and the classical 2D system is then completed by using: 
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where y is the vortex fugacity. In the 2D XY model (1) 
(with Jc = 0) one has a single energy scale given by Jab, 
so that the vortex-core energy n is given by:[l, 2] 



^J'XY = TT Jah ln(2\/2p7) ~ LGirJab 



(6) 



where 7 is the Euler's constant. However, /x depends in 
general on the details of the microscopic superconducting 
model under consideration, so it will be taken as a free 
parameter in the following, while the value (6) will be 
used just for the sake of comparison with the XY model 
(1). It is worth noting that the limitations of the XY 
model as an effective phase-only model have been pointed 
out in Ref. [18], as far as the role of the phase- interaction 
terms beyond Gaussian level are concerned. The effect of 
the interlayer coupling of Eq. (1) can be incorporated 
in the sine-Gordon model (3) as an interchain hopping 
term, so that the full Hamiltonian becomes: 
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where m is the chain (layer) index and gj^ = ttJc/T. We 
derived the perturbative RG equations for the couplings 
of the model (7) by means of the operator product ex- 
pansion, in close analogy with the analysis of Ref. [21]. 
Under RG flow an additional coupling g± between the 
phase in neighboring chains is generated: 
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') + -^{dx(l>7n){dx4'm') 
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The superfluid coupling Ks is deflned, as usual, [17] as the 

second-order derivative of the free energy with respect an 
infinitesimal twist 6 of the phase, dx9m dxdm — 5. The 
interlayer term (8) contributes bs K ^ K{1 — ng±) to 
the current-current coefficient, where n = 2 is the number 
of nearest-neighbors chains (layers). Thus, the in-plane 
stiffness J, is defined as: 
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The full set of RG equations for the couplings 
K,Ks,gu,gjc reads: 
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with e = log{a/ao), where ao,a are the original and 
RG rescaled lattice spacing, respectively. Observe that 



3 



for gj^ = the first two equations reduce to the stan- 
dard ones of the KT transition, [1] with a fixed point at 
K — 2,gu ^ 0, and Kg coincides with K. Thus, one sees 
that at K > 2 the f/„ couphng is irrelevant, the quan- 
tum ID system is a LL and the vortex-antivortex pairs 
are bound in the classical 2D system. At K < 2 the g„ 
term is relevant, the field is locked in a minimum of 
the cos(2^) potential and the ID system is an insulator. 
In the classical case this corresponds proliferation of vor- 
texes (large vortex fugacity) in the metallic phase. The 
physical superfluid stiffness is given[2] by the asymptotic 
value of the running couphng Ks{£) = K{£). Thus, Jg 
is finite below Txt, since K{£) flows to a finite value (in 
particular K{oo) = 2 at T^t, in accordance to Eq. (2)), 
and it goes to zero above the transition, since K{£) scales 
to zero. The KT temperature is defined by the highest 
temperature where K{oo) = 2, and it is given (at small 
gu) by K{T) ~ 2 = 2g„(T), which yields Tkt ^ Jab/2. 
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FIG. 1: (Color online) RG flow of the couplings K(i) and 
Ks(€) at various temperatures for = [i,xy and r\ — 10~*. 
Inset: critical temperature Tc as a function of r\ for a bare 
stiffness Jo(r) = Jab(t ~ T/AJab) (see text). 

As an initial value gj^ ^ is considered, the interchain 
(interlayer) coupling increases under RG, [14-16] leading 
to larger values of the LL parameter K{t} and stabiliz- 
ing the metallic ID phase. However, when the initial gu 
coupling is sufficiently large the second term in the r.h.s 
of Eq. (10) dominates and K{t) goes to zero, leading to 
the insulating ID phase. In the classical 2D analogous 
the effects of gj^ are easily readable trough the behav- 
ior of Ks, which is controlled by the gu coupling alone. 
Whenever K{t} scales to large values the gu coupling is 
irrelevant and Kg flows to a constant, see Fig. 1. This 
effect guarantees the persistence of the superfluid phase 
in a range of temperature above Tkt- Indeed, the initial 
decrease of Ks{i) is cut off at a finite length scale by the 
interlayer coupling, which brings again K{£) to large val- 
ues and guii) to zero, giving a finite asymptotic value of 
Ks{i). As the temperature increases further and the gu 
term dominates, both K and Kg scale to zero, the scal- 
ing dimension of gj^ becomes negative and one observes 



a "layers decoupling" above Tc. [15, 16] 
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FIG. 2: (Color online) Temperature dependence of Js{T) in 
2D (lines) and in the layered 3D case (symbols). Here Jo{T) = 
Jab{l — T /AJab)- The Tkt is identified by the intersection 
between Js{T) and the straight line 2T/it. The results for 
M ^ Mxy show a rapid downturn of Ja{T) at Td — Tkt- As 
increases Td increases as well, so that at Tkt no effect is 
observed in Js{T) reminiscent of the jump present in 2D. 

The critical temperature T^ is defined by the vanish- 
ing of the Ks{£ oo). Alternatively, to account for the 
perturbative character of the RG equations, one can com- 
pute Js by stopping the RG flow at the scale £* where 
gj^ is of order one. [23] The two deflnitions are equivalent, 
and lead to the estimate of the critical temperature Tc 
reported in the inset of Fig. 1. For the sake of complete- 
ness, we also added a temperature dependence of the 
bare couplings, using Jo(T) = Jafc(l — T/AJab), as due 
to long-wavelength phase fluctuations in the XY model 
(1),[17, 18] and we keep the ratios rj = Jd Jab and [ij Jah 
fixed. For [i — [Ixy the calculated values of T^j Jab show 
a remarkable quantitative agreement with Monte Carlo 
simulations on the anisotropic XY model. [17] The fact 
that larger values of lead to a larger critical tempera- 
ture has a direct counterpart on the temperature depen- 
dence of the superfluid stiffness Js (T) , as we show in Fig. 
2. As one can see, when Jc = we recover the standard 
jump of Js at Tkt, which is easily identified as the tem- 
perature where the curve Js{T) intersects the line 2T/7r, 
according to Eqs. (2) and (9) above. As n increases the 
gu coupling decreases, and the renormalization of Js{T) 
with respect to Jo{T) below Tkt becomes negligible. As 
soon as a finite interlayer coupling is switched on, the 
jump of Js{T) at Tkt disappears and it is replaced by 
a rapid bending of JsiT) at some temperature Td- How- 
ever, while for /i < fixY Td coincides essentially with 
Tkt, for a larger vortex-core energy Td rapidly increases 
and approaches the temperature Tc estimated above. 

These results offer a possible interpretation of the 
experiments in underdoped YBC0.[12]], where rj ^ 
10~^[20], but the measured Js{T) goes smoothly across 
the Tkt estimated from Eq. (2). We calculated Js{T) as 
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FIG. 3: (Color online) Correlation length above Tc for dif- 
ferent values of rj and ^. Each curve stops at Tc, which in- 
creases as rj increases. We take Jo{T) = Jab{^- — T /Tmf), with 
Jab = 300 K and Tmf = 100 K, as appropriate for Bi221.[8] 
The dashed line is ^ ~ g/3M(T)/2^ expected to hold far above 
Tc- The solid lines are fit with the KT functional form (see 
text), with Td = Tkt,c = 0.25, fe — 0.9 for /i = hxy and 
Ta = 86.8 K, c = 0.4,6 = 1.21 for ^ = SfixY- 



anisotropic layered SC. Using a RG approach, we showed 
that a finite interlayer coupling can shift the temperature 
scale Td of vortex unbinding away from the KT temper- 
ature Tkt of the pure 2D case. Indeed, Td is essentially 
controlled by the vortex-core energy, and it coincides with 
Tkt only when /i < fJ-xY, as within the standard XY- 
model (1). When applied to cuprates, our findings sug- 
gest that in these systems /i is definitively larger than ex- 
pected in the XY model, even though still of order of the 
in-plane stiffness. The consequences are twofold. First, 
the lack of any signature of KT behavior in Jg (T) at Tkt 
does not rule out the possibility that phase-fluctuation ef- 
fects play a role in these systems. Second, fi > ^xy is not 
inconsistent with microscopic theories which associated fj, 
to the energy scale of the superfluid stiffness instead of 
that of the superconducting gap, which would be far too 
large compared to Jg in underdoped samples. [3] 

We ackonowledge useful discussion with S. Caprara, 
M. Cazalilla and M. Grilli. This work was supported in 
part by the Swiss NSF under MaNEP and division II. 



done in Fig. 2 (taking into account also the measured lin- 
ear depiction at low T) . Using a large vortex-core energy, 
i.e. /i — GfixY, we found that Js{T) shows no signature 
of a rapid downturn at Tkt, and goes to zero near to the 
measured Tc- Observe that we did not consider the effects 
of disorder, which can also smear out the KT transition, 
as measurements in thin films[4, 25] could suggest. 

A similar separation between Tkt and Td is observed 
in the behavior of the correlation length above Tc. Since 
the quantity experimentally accessible is the vortex den- 
sity ny, given at the RG scale £ = ln(a/ao) by nv{i) = 
g-/3M(^)/Q,2 _ g„(i!)/(27ra^), we define the correlation 
length ^ as = nv{is), where Us — aoe^= is the 
length scale where Ks{is) vanishes above Tc- The be- 
havior of ^ (T) for different values of /i and rj is reported 
in Fig. 3, using parameter values appropriate for Bi2212 
compounds. Far above Tc ~ ao, so that ^ scales as 
flo/ ■y/(?„(T)/27r, as shown by the dashed line. As T ap- 
proaches Tc is increases and shows the exponential in- 
crease reminiscent of the KT behavior in 2D.[1] However, 
while in 2D ^ diverges at Tkt, a finite Jc cuts off at Tc 
the increasing of ^, since below Tc gu becomes irrele- 
vant and Kg flows to a finite value. Nonetheless, the 
behavior of ^ above Tc is still reminiscent of the KT be- 
havior, S^kt ~ floc exp(6/ \/T jTd — 1), with c, h of order 
one, provided that Tkt is replaced by a proper scale Td 
slightly smaller than Tc- Once again, while for \x — [Ixy 
Td ~ Tkt, as the vortex-core energy increases Td be- 
comes significantly larger than Tkt and approaches Tc. 
This behavior is consistent with recent experiments in 
Bi2212 compounds. [7] 

In summary, we analyzed the phase-fluctuations 
contribution to the 2D KT-3D crossover in strongly 
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